Static Solutions with Spherical Symmetry in f(T) Theories 



Tower Wang* 

Department of Physics, 
East China Normal University, 
Shanghai 200241, China 



(Dated: January 12, 2013) 



Abstract 

The spherically symmetric static solutions are searched for in some f(T) models of gravity 
theory with a Maxwell term. To do this, we demonstrate that reconstructing the Lagrangian of 
f(T) theories is sensitive to the choice of frame, and then we introduce a particular frame based on 
the conformally Cartesian coordinates. In this particular frame, the existence conditions of various 
solutions are presented. Our results imply that only a limited class of f(T) models can be solved 
in this frame. For more general models, the search for spherically symmetric static solutions is still 
an open and challenging problem, hopefully solvable in other frames. 
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I. INTRODUCTION AND PRELIMINARIES 



As a great triumph of science in the last century, Einstein's theory of general relativity 
with a cosmological constant has been established by many observations and experiments 
hitherto. In spite of this fact, people have been enthusiastically trying to modify or extend 
this theory. Cosmologically the motivation to do it is running our Universe with less dark 
matter Ej2| or without the cosmological constant or driving inflation without a scalar 

field [(1,17] or with richer phenomena jsj. 

A few years ago, f(T) theories of gravity were proposed as an alternative of the cosmo- 
logical constant to explain the accelerated expansion of the late Universe Similar to 
the f(R) theories, f(T) theories deviate from Einstein gravity by a function f(T) in the 
Lagrangian, where T is the so-called torsion scalar. This class of theories has received a lot 



of attention recently [lCH30j , but earlier examples could be traced to (3l|, |32| . In order to 
better study /(T) gravity theories, we will seek for exact solutions of the field equations, 
akin to the work |33l-l37j done in f(R) models. Especially we are interested in spherically 
symmetric solutions. The widely studied Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric is such a solution obviously, but we will pay attention to static solutions instead. 

In this paper, the notation of indices is as follows: Greek indices /i, u, ... and Latin indices 
from the beginning of the alphabet a, b, ... run over 0, 1,2,3, while Latin indices from the 
middle of the alphabet i, j, ... run from 1 to 3. In our notations, the Kronecker delta is defined 
with the same normalization irrespective of the position of indices, e.g. 5qq = 5 00 — 5® — 1. 

The dynamic variables in f(T) theories are the tetrad (vierbein) fields e a . In a coordinate 
basis, the components of dual tetrad fields e^dx^ are related to the familiar metric tensor 
g^ through 



1 



where rj ab = diag(l, -1, -1, -1) 
to construct the torsion tensor 



Starting from the tetrad fields, one may follow |9|, |1 



the contorsion tensor 



and a useful tensor 



(2) 
(3) 
(4) 



Then in terms of the torsion scalar 



and the notation e = det(e^) 



I = J d 4 xe 
which leads to the equations of motion 



g, the action of f(T) gravity theories is given by 
1 



167rG 



(T + /) + £, 



(5) 



(6) 



;i + f, T )[e-%(eSr) ~ e^ A S/1 + IttS^T - -<(T + f) = ^Ge'T/ (7) 
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with SJ 1 " = e^S ^ v . If the function f(T) is replaced by a constant, the theory is equivalent 
to Einstein's theory with a cosmological constant (24j . 
Throughout this paper, the Maxwell term 

C m = ~V F ( 8 ) 
will be considered. It gives the energy-momentum tensor 

77 = F pX F» x - \5»F X »F X ». (9) 
We will need the Ricci scalar R = g^R^, taking the convention of Ricci tensor 

Rfiv = d\T^ u — d u v^ x + r Ap r£„ — r^ p r^ A (10) 

with the Levi-Civita connection 

rjL = \9 Xp {^g pv + d u g, p - d p9pM ). (11) 

In the rest of this paper, after choosing the frame in section [Til we wm search for the 
spherically symmetric static solutions in different situations in sections \IV\ IVl Then we 
will recap the solutions and their constraints on Lagrangian in section IVIl In section IVIII 
some open issues will be mentioned. 



II. THE CHOICE OF FRAME 

The first difficulty we encountered is the choice of frame. As dynamic variables in f(T) 
theories, tetrad fields are sensitive to the frame. More importantly, the torsion scalar is also 
frame-sensitive, very much unlike the Ricci scalar. This trouble is attributed to the lack of 



local Lorentz invariance in f(T) theories 24]. Consequently, even for the same metric and 



the same coordinate basis, different frames result in different forms of equations of motion. 

To see this point clearly, we take the FLRW metric for example. Corresponding to the 
coordinates 

ds 2 = dt 2 - a 2 (t)5 ij dx i dx j , (12) 
the following form of tetrads are widely used in the literature 0, [Io[ , 

= ajdz", e^dx^ = otf'dz". (13) 

It give the torsion scalar T = —GH 2 and the Ricci scalar R = —6(2H 2 + H). Here the 
Hubble parameter H = a/a, and the dot overhead implies the derivative with respect to t. 
Then equations of motion (j7j) lead to the generalized Friedmann equations 

/ + QH 2 + YlR 2 f? = lQnGp, 
f + QH 2 + AH + 4(3# 2 + H)f, T -48H 2 Hf :TT = -lQnGp. (14) 

Expression ([TBI is an obvious form of tetrads yielding metric (|12|) . But it is not the unique 
choice. For instance, we can rotate it to a different frame, and write down the following 
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form of tetrads: 



e°cb? 

fj, 

eldx 1 " 



eldx" 



(it, 

-(xdx + yrfy + zciz), 



V^ 2 + ?/ : 



\fx r Ty- 



■.{xdx + ydy) - 
— xdy), 



x 2 + ?/ 2 



dz, 



(15) 



+ y 2 + z 2 ) 1 / 2 . It is easy to check that the new form of tetrads can also 



in which r — ~ 
produce metric (fT2l) and the local rotation matrix 



/ 



v 

r 

yz 



ry/ x 2 +y 2 ry/ x 2 +y 2 

_y x 

a/ x 2 +y 2 



o 



(16) 



/ 



transforms ({13]) to (|T5|) . Keep in mind that the rotation group is a subgroup of the proper 
Lorentz transformation, so frames (|T3|) and (|T5|) are related by a local Lorentz transforma- 
tion. In accordance with tetrads ( !T5|) . the torsion scalar becomes T = 2a~ 2 r~ 2 — 6H 2 but 
the Ricci scalar remains R = —6(2H 2 + H). Now equations of motion ([7]) reduce to 



f + 6H 2 + l2H 2 f >T 



f + 6H 2 + AH + 4(3# 2 + #)/ )T 



a?r 2 



f. 



T 



f. 



TT 



a?r 2 



16nGp, 
0, 

— lQnGp. 



(17) 



Comparing (|T4l) and (|T7|) . we can see the equations of motion, as well as the reconstruction 
of f(T), is quite sensitive to the choice of frame. Especially, the second equation of (JT7J) 
requires / = AT — 2A. After rescaling the Newtonian constant, the resulted theory is simply 
equivalent to Einstein's theory with a cosmo logical constant. But the commonly studied 
equations f[T4"l) . which are derived in frame f|T3|) . allow much more general forms of /(T) 
to survive. This property does not rely on the choice of coordinate system 1 . Through the 



coordinate transformation x 
traded to 



r sin 9 cos 0, y — r sin 9 sin < 



e° u dx" 



dt, e^dx^ = adr, 



e 2 u dx" 



ard9, 



z = v cos 9, expression ffT5]) is 
ar sin 9d<p (18) 



eldx* 



in a spherical coordinate system 

ds 2 = dt 2 - a 2 (t)(dr 2 + r 2 dn 2 ) (19) 

with dfl 2 = <i# 2 + sin 2 9dcj) 2 . That is to say, (|T5|) and ({TBI) describe the same frame in different 
coordinate bases. 



The author is indebted to Rong-Xin Miao and the referee for their valuable comments on this point. 
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The above example sharpens the problem of what frame we should adopt to look for new 
solutions. It is interesting to notice that the spatial coordinates in metric (Tl2|) are Cartesian 
up to a conformal factor. This coordinate system is relatively more convenient for us to 
perform the Lorentz boost and the rotation of spatial axes locally. Fortunately, there is a 
close cousin of (1121) . namely the generalized Gullstrand-Painleve metric 

= + (20) 

where a and (3 are functions of the radial coordinate r = (5ijX z x^) 1 ^ 2 for static solutions. 
In this coordinate system, every x^ fully spans (— oo, +oo), and the spatial coordinates are 



conformally Cartesian [38|, |39| . Through a coordinate transformation 



metric (120]) can be rewritten in the more familiar form 

y/adt + -j J - ^pdtt 2 

in spherical coordinates. 

From the metric in Gullstrand-Painleve coordinates, we directly read out one possible 
form of the tetrad fields 

e°dx" = dt, eidx" = ^-dt + \dx\ (23) 

This is the closest analogue of ffl3|) for the solution we are interested in. It is also the starting 
point of our investigation. The impact of such a choice of frame is twofold. On the one hand, 
it simplifies our calculation significantly, otherwise the Lorentz transformation will induce 
six more variables and then the equations of motion would be too cumbersome to tackle. 
On the other hand, it restricts the validity of our solutions to a limited class of f(T) models, 
while the search for spherically symmetric static solutions in more general models remains 
a challenging task. 



III. HUNTING FOR SPHERICALLY SYMMETRIC STATIC SOLUTIONS 

In the Maxwell term (JS}, F^ u = d^A v — d„A^ for an electromagnetic field. If the field is 
static and spherically symmetric, we can make use of U(l) gauge invariance to write it in 
the form A^dx^ = ^f{r)dt. The electric charge of our solution can be obtained by evaluating 
the integral 

1 

in the limit r — > oo. 



-99 9 



tt " rr F tr d6d V 



r 1,; 



(24) 
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Starting from tetrads of the form (1231) . we find the torsion scalar 

T = ~W ~ r P,r) 2 + rfta r (ft - rft >r ) - r 2 ftf r }. (25) 

The equations of motion are complicated, whose tedious form will not be shown here. We 
succeeded in casting them in simple forms under different situations. One situation is that 
with a = 0. This will be studied in section |TVJ The other situation is to be investigated in 
section [V] where a^O. 



IV. a = 0, ft ^ 



For physical solutions, ft should be non-zero, but a = is still allowed. In this situation, 
the torsion scalar simplifies as T = 2 ft' I, and the equations of motion reduce to 



rft 



A. 



;i + h 



T_2ftft, 
2 r 



;i + fa 



(ar - (i + f, T ) -(r + f) + MGftW, 



rft 



A. 



2_,2 



+ 8ttG73 2 7 



2 

7,1 



0. 

0. 
0. 
0. 



Equation ( |26l) can be integrated, resulting in the relation 



A 



+ /,: 



(26) 

(27) 
(28) 

(29) 
(30) 



Here Q and are constants, whose physical significance will be clear later. We will assume 
Too > without loss of generality. Inserting this relation into (128"]). we get 

SpF- 2 ) +rV ' = ° (31) 

because /3 7^ 0. This equation can be combined with 
differential equation of 7, 

5! , 3 

2 7,rr + ^7,r 



0. 



which leads to 



to give a second-order ordinary 

(32) 

rooV r + r 

with a nonnegative constant r . We have chosen a proper signature to ensure that Q is the 
electric charge given by f )24|) . Now the 7 2 term can be eliminated in equation (|3T|) and the 
equation becomes 

Q 2 frft tr r 2 + 2r\ 



7,! 







ft r 2 + r 2 J °' ^ 

When Q/tqo 7^ 0, this is a first-order equation of ft. We will continue our discussion in two 
situations: Q/r^ = in subsection IIV Al and Q/r^ 7^ in subsection IIV Bl 
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A. Q/roo = 

In the situation Q/ r oo = 0, equation fl33l) indicates 7 = up to a constant eliminable by 
U(l) gauge transformation. At the same time, the right hand side of equation (1301) vanishes. 
To this equation there are two solutions. 

One solution is a constant /3, and we fix it to be (3 = 1 by rescaling the spatial coordinates 
x l . This solution, summarized as 



a = 0, 

(3 = 1, 

7=0, (35) 

describes a Minkowski spacetime. One may check the torsion scalar T = and the Ricci 
scalar R = 0. It is easy to also check that the constant solution satisfies fl2"B"l) . f[2"Sj) and 
(T25j) automatically. But the solution is incompatible with equation (12?]) unless /(0) = 0. 
That does not mean the gravity theory must recede to the Einstein's theory to guarantee 
the existence of global Minkowski solution. In fact, it provides a sufficient but not necessary 
condition: the global Minkowski solution exists if function f(T) vanishes at point T = 
numerically. If this condition is violated, the Minkowski solution might still exist in frames 
other than (1231). 

The other solution to (1301) is 

1 + f,r = 0. (36) 

Then equation (|27|) tells us 

T + f = 0. (37) 

If T is not a constant, one should solve equations (1361) and (137]) analytically as differential 
equations. The only consistent solution to them is / = — T analytically. However, this 
solution is unphysical and should be ruled out, otherwise the theory of gravity is ill-defined 
by a null Lagrangian. If T is a constant, it will be enough to require these equations to hold 
numerically, not analytically. In other words, we do not have to solve them as differential 
equations. Remembering that T = 2/3 2 r , a constant T can be achieved by 

a = 0, 
r 

f3 = — + k, 
'00 

7=0. (38) 

Here the constant k = or ±1 after rescaling x l . From ( l38l) we work out the torsion scalar 
T = and the Ricci scalar 

a = -4 - - (») 

Clearly there is a curvature singularity at r = unless is infinite or k = 0. So let us 
consider different choices of parameters. 

1. Too is infinite and k = 0. This solution gives (3 = and an ill-behaved metric. 

2. Too is infinite and = ±1. We identify this solution as the Minkowski solution fl35|) . 
but emphasize that equation fl36l) is not a necessary condition for the existence of 
Minkowski solution, as shown at the beginning of this subsection. 
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3. Too is finite and k = 0. In terms of a new coordinate z = r M ln(r/roo), the metric from 
this solution takes the form 

ds 2 = dP- dz 2 - r^dQ 2 . (40) 

This metric describes a spacetime of topology R 1 x R 1 x S 2 , where the S 2 sphere is 
of constant radius r^. 

4. Too is finite and k = ±1. It is easy to see the curvature singularity at r = is a naked 
singularity in these solution. 



B. Q/r-oo ^ 

In situations with Q/r^ ^ 0, equation (13"4"1) is solved by 



,2 



= 7t=^' ( 41 ) 



(42) 



which immediately gives 

T _ 2r 2 (r 2 + 2r 2 ) 2 

r 2 o ( r 2 + r ,2^3 ' 

One may use (124")) . (1331) and (j4"Tl) to prove that indeed Q is the electric charge. 
Substituting equations fl33]), (1M]1 . (EE!]) and (122"]) into ((27} and ([30]), we obtain 

r = 87rGQ 2 r 2 (r 2 + 4r 2 ) 
r^(r 2 + r 2 ) 2 
SvrGQV + r 2 ) 
1 + / ' T ~ ^(r 2 + 2r 2 ) • (43) 

Making use of fl4"2"]) . one may check the consistency condition (T + f) <r = (1 + f,r)T^ by 
straightforward calculations. 

In order to work out more details, we should consider two possibilities. 

1. ro = 0. For such a choice, we get the solution 

a = 0, 
r 

= — , 

7 = -^ln(^V (44) 



In the notation z = ln(r/roo), it gives metric (1401) and 7 = —Qz/r 2 ^. Furthermore, 
equations ( 1431) become 

8irGQ 2 



T + f 



r 4 

00 



1 + It = —5—, (45) 

00 

which should hold numerically for T = 2/r 2 ^. Extrapolated to the limit Q = 0, 
this solution reproduces the regular R 1 x R 1 x S 2 solution presented in the previous 
subsection. 
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2. tq > 0. Then as a cubic equation of the variable r 2 /r^, equation (jJ2J) has one real root 
r 2 jr 2 = u(T) and thus 



T + f 



8nGQ 2 u(u + 4:) 



1 + / ^^±il. (46) 
r 2 (w + 2) 



The metric is decided by the solution 

a = 0, 

P = 



^ooV /r2 + r 



7 = — arcsinh ( — J — ^-arcsinh ( — ) . (47) 

^oo V r 0/ ^oo V r 0/ 



In terms of coordinate 



arcsinh ( — ) — arcsinh 




(48) 



we write it as 



ds 2 = dt 2 - dz 2 - r 2 dQ 2 , (49) 
which describes a spacetime of topology R 1 x R 1 x S 2 . The radius of S 2 sphere is 



r 2 



1 + -f. (50) 



In the region r G (0, oo), both z and f are strict monotonic functions of r, so the 
radius f shrinks monotonically as z varies from — oo to +oo, as illustrated in figure [TJ 
Extrapolating this solution to the limit r = 0, we recover equations (j40p . f|44|) . ( I45p 
as well as T = 2/r 2 ^ and z = ln(r/roo). Note the extrapolation from f T46|) to f )45|) is 
a little tricky: the special metric ( 14"0"|) exists as long as equations 045p hold numerically 
at the point T = 2/r 2 ^, but f j46l) should stand analytically to guarantee the existence 
of the general metric ( 1491) . 
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Figure 1: The radius f of S (|50j) decreases monotonically as coordinate z ([45]) increases. 
V. a/3 / 

When a/3 ^ 0, the equations of motion can be arranged in the form 

(1+/ ,)(r^ = o, (5D 

f, TT T r = 0, (52) 

^2T - ^) (1 + /,t) ~(T + f) + 87rG/3 2 7 2 r = 0, (53) 



2a - (3^ + >r - a >r p - 2/3 ^ - a iTr P' 



(1 + f >T ) - 167rG/3 2 7 2 r = 0, (54) 

= 0. (55) 



The last equation is the non- vanishing component of Maxwell equations. 

One special case is /t = —1. Keep in mind that (3^0, then equation ( 154"]) indicates 
7 = up to i7(l) gauge transformation. 

If f t x 7^ — 1, equations (1ST]) is solved by 



" = UJ • (56) 

where we have chosen a positive signature for without loss of generality. According to 
equation ( 155]) we can write down the relation r 2 7 r / /3 = —Q with a constant electric charge 
Q, as demonstrated by ( I24p . This relation implies 

7 = -^lnf-l (57) 



for p = 1 or 



(1 — p)r \r 

for p ^ 1 up to an eliminable constant. 
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We will study the case with f t T = — 1 to some extent in subsection IV Al The case with 
f t T 7^—1 will be exhaustively investigated, for p = 1 in subsection IV B\ and for p ^ 1 in 
subsections |VC] (T r ^ 0) and EE (T = const). 



A. f, T = -1 

In this case, the equations of motion are further simplified, 

1 + f,T = 0, 
f,T T T r = 0, 

T + f = 0. (59) 

If T r 7^ 0, equations ( 159"]) have only one solution / = — T, resulting in an unphysical null 
Lagrangian. So we turn to consider T = const. 

The expression of T is given by f l25|) . It is difficult to exhaust all of the solutions to 
T r — 0, but we found a particular solution 

r 

P = — 

'"oo 

7 = 0, 

1 + /,T = 0, 

T + / = 0, (60) 



leaving a unconstrained. The corresponding metric is 



,j. 2 



ds 2 = (1 - a)dt 2 - (1 _^ )r2 ^ 2 - r L^ 2 - (61) 

This metric generalizes the R 1 x R 1 x S 2 solution (|40|) of constant radius, but the constraints 
on Lagrangian are the same. For this solution, the torsion scalar T = and the Ricci 

scalar R = — (r 2 a, rr + ra )T + 2)/r^ . 
Another particular solution to fl59l is 

Ar 2 2GM 

« = -ir + ' 

3 r 

P = 1, 

7 = 0, 

1 + /,t = 0, 

T + / = (62) 

with T = — 2A. Here A may be thought as the "cosmological constant" though this model 
does not include the Einstein's gravity theory. Later in subsection IV Dl we will see such 
a solution actually exist for more general models with any value of f t T, which naturally 
incorporate the Einstein's theory. Please refer to subsection IV Dl for physical interpretations 
of this solution. 
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B. f iT ^ -1, p = 1 



In the previous subsection, we have studied some solutions with = —1. From now 
on, we will consider the case of f t T ^ — 1- In this case, one special example is p — 1. Now 
the solution of is given by ( loT)|) with p — 1, while the solution of 7 is ( 15"T|) . Then one may 
directly check that T = 2/r 2 ^ using equation f[2"5"j) . Subsequently, most equations of motion 
are automatically satisfied except for ( lol?|) and (j5"3|) . which demand 



8nGQ 2 



(63) 



and 



1 



8nGQ 2 



^wj^U^U +c (64) 

respectively. In the above solution, g and C are constants, and we have taken into account 
the fact that is independent of r. It is easy to show that the corresponding geometry 
contains an S 2 sphere of constant radius r^. More interestingly, solution (jUJ) and constraints 
(jlSJ) can be regenerated from this solution by requiring a = 0. 

In what follows, restricting to the case 7^ — 1, we will explore the possibility that 
p 7^ 1. We will put our results in two subsections, according to whether the torsion scalar is 
a constant. 



In this case T is not a constant, so the solution to equation flo^]) is T+f = (T — 2A)/A. In 
order to give the correct value of Newtonian constant, A = 1 should be imposed on physically 
viable models. The viable model is simply the Einstein's gravity theory with a cosmological 
constant A. Although the model is not new, in this subsection we will briefly cover it for 
completeness. 

Since f(T) is a constant, from equations (153!) and (|54|) we obtain 
T = -2A + -77 ( — " ' 



r 2 \ r oo7 r A vr-oo, 

p(p-2) | Cr 2 /r^yp | 2GM /rV AttGQ 2 f r \ 2p 

(p — l) 2 (p — l) 2 V r / (p — l) 2 r \ r oo/ (p — l) 2 r 2 \rocJ 

Here M and C are constants of integration, while G is the Newtonian constant. Comparing 
the solution with (125]) . we can identify C = A/3. 

With the redefinition f = r^r 1 "^, this solution can be reformed as 

Af 2 2GM A-nGQ 2 



1-a = T77 1 ^ + 



(p — l) 2 \ 3 r r 2 

/3" 2 r 2 = f 2 , 

(3- 2 dr 2 = - 1 , 9 rff 2 , 



7 



(p-1) 2 



(1 — p)f' 

/ = -2A. (66) 
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After the time coordinate is rescaled as t — > (1 — p)t, solution fl66|) recovers the familiar 
Reissner-Nordstrom metric, which describes a charged black hole in de Sitter (A > 0) or 
anti-de Sitter (A < 0) or flat (A = 0) spacetime. Other than the coordinate transformation 
in this way, we can also recover the metric directly from (15B1 . (158]) and fl65|) by setting p = 0. 
In this solution, the integration constants M and Q are mass and electric charge of the black 
hole respectively, while A is the cosmo logical constant. 



D. f )T ^ -1, p 1, T r = 

Since p ^ 1, the expressions of (3 and 7 are given by equations (15B1 and (1581 . We insert 
them into equation fl53|) and expand the resulted equation in powers of r, 



^ T . 1+T .% (l+M (i\\^f(i\''^ (67) 

H \ r oo/ V r oo/ 

The first two terms on the left hand side are constants, because / and f t T are constants 
when T takes a constant value. Remember also that 7^ — 1, P 7^ 1, so this equation 
dictates p = 0, Q = and 

2T/> - / + T = 0. (68) 
Therefore, the consistent solution is restricted to 



P = l, 7 = 0. (69) 

Let us denote T = — 2 A as a constant. Then equation (1251) can be solved as a differential 
equation of a by 

Ar 2 2GM 

a = — + . (70) 

6 r 

One may check that equations ( 1521 and (J5i|) are satisfied automatically. 

Extrapolating ( l68l) . ( l69l) and ( 1701 to the limit f t T = — 1, we reobtain equations (|62l . thus 
equation (16"8"|) holds more generally as an existence condition of this solution. The solution 
can be interpreted as a Schwarzschild black hole in de Sitter (A > 0) or anti-de Sitter 
(A < 0) or flat (A = 0) spacetime. The integration constant M is related to the black hole 
mass, while A is a constant determined by the Lagrangian. Although A does not necessarily 
arise from a constant term in the Lagrangian, we will still call it cosmological constant for 
simplicity. It is straightforward to check that the Ricci scalar R = — 4A. 

We observe that torsion scalar T and equation (l6"8l depend only on "cosmological con- 
stant" A, so the mass parameter M does not enter into Lagrangian parameters. This is 
consistent with our physical interpretations of the integration constants. Equation (168]) puts 
a constraint on Lagrangian of viable models that admit such a solution. Since T is a con- 
stant, it is enough to require this equation to hold numerically, not analytically. In other 
words, we do not have to solve it as a differential equation. 

In subsection IIV A\ we have found out a solution describing the Minkowski spacetime. In 
frame (|23p . the Minkowski solution exists in f(T) theories satisfying /(0) = numerically. 
Actually, as can be confirmed directly, this solution and its constraint on Lagrangian can be 
obtained as a limit M = 0, A = of the Schwarzschild solution in this subsection. 
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VI. RECONSTRUCTION OF LAGRANGIAN 



The results in sections IIVI and |V] look scattering. In the current section let us review 
them briefly by collecting the existence conditions of some solutions. In other words, we 
will review the conditions f(T) should meet if certain solutions exist in frame (1251) . As will 
be shown, some of the results can be unified in the same form, although they were derived 
under different assumptions. As a convention, when we say an existence condition should 
hold "analytically", we mean the equations should stand for any value of T, solved as a 
differential equation analytically. Otherwise, the word "numerically" is used to mean that 
the equations are only to hold for a specified value of T. We should stress that all of the 
existence conditions here apply only to the frame (]23"1) . But we do not claim anything about 
other frames. Particularly, when the conditions are violated, these solutions may appear in 
a different frame. 



A. Reissner-Nordstr6m-(anti-)de Sitter Solution 



We have investigated this solution in subsection IV CI Using the notations of (1231) . this 
solution is 

pip -2) Ar 2 fr^p 2GM ( r \ p AttGQ 2 fr^ 2p 

« = ~, 7777 + 777 7T77 + " 



(p — l) 2 3(p — l) 2 V r / (p — l) 2 r \TooJ (p — l) 2 r 2 \r c 

P = 











Q 


/ r 


(1 — p)r 


\roo 


Solution 


(I7TD 



7 = 77^: - (71) 



/ = -2A (72) 



analytically, which is equivalent to the Einstein's theory with a cosmological constant A. We 
achieved at this condition in subsection IV CI under the assumption ^ — 1- Nevertheless, 
by combining the fact that T r ^ for this solution and the discussions at the beginning of 
subsection IV A| we can see the condition is robust even without any assumption on f\x- 



B. Schwarzschild-(anti-)de Sitter Solution 



The solution is 



Ar 2 2GM 

a = — + , 

6 r 

P = 1, 

7 = (73) 



in frame (J23J). The existence condition of solution (173]) is 

2T/>-/ + T = (74) 
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numerically with a constant T = — 2A. We can also consider f|74l) as the existence condition 
of (anti-)de Sitter or Minkowski solution in that frame, depending on the value of A. This 
solution has been studied in subsection IV Dl which stands on a general ground. The special 
case with fp = — 1 was explored in subsection IV Al while the limit of M = 0, A = was 
investigated in subsection IIV Al 



C. Charged Solution with a Sphere of Constant Radius 

This solution is 

1 - 



a. 







7 



8nGQ 2 



1 + /, 



T 



+ In 



+ c, 



in frame (123]) . Its existence condition is 

T + f 



8ttGQ 2 



(75) 



(76) 



numerically, where T = 2/r 2 ^. The above solution was obtained in subsection IV B\ but the 
result can be extrapolated to the special case a = 0, namely the R 1 x R 1 x S 2 solution in 
subsection IIV Bl In that limit, there is an additional existence condition 



1 + /. 



8nGQ 2 



T 



(77) 



This condition can be successfully recovered from the requirement a = 0. Further taking 
limit Q/Vqo = 0, we find the neutral R 1 x R 1 x S 2 solution discussed in subsection IIV Al 
The neutral solution can be extended in another way, as to be elucidated in the coming 
subsection. 



D. Neutral Solution with a Sphere of Constant Radius 



In frame (1231) . the solution is 

r 

fi = — 

'"oo 

7=0, (78) 

leaving a unconstrained. This solution exists in frame (T23"]) if the Lagrangian satisfies con- 
ditions 

1 + f,T = 0, 

T + f = (79) 

numerically at T = 2/r 2 ^. We have investigated this solution in subsections IIV Al (a = 0) 
and IV Al (a / 0). Both solutions (JT5j) and (1781 have a product space S 2 of constant radius 

Too- 
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E. Solution with a Sphere of Variable Radius 



As has been discussed in subsection HVBt the solution 

a = 0, 

= 



7 = — arcsinh ( — J — — arcsinh ( — J . (80) 
V r o/ r oo \r J 

has a product space S 2 of variable radius (1501) . The existence condition of this solution in 
frame (I23p can be cast in the form 

T W M + 4) 

r^(w + l) 2 

which is expected to hold analytically. In condition function w(T) is implicitly given 
by the real root of cubic equation 

2n(^ + 2) 2 

r^u + l)*" 1 ) 

Note the second equation of f|4"6l) can be derived from the first one and equation fl82l . 



F. Example of Lagrangian Reconstruction 



The existence conditions summarized above are useful. On the one hand, given the 
Lagrangian of /(T) model, we can quickly judge the existence of these solutions in frame 
(123]) . On the other hand, we can use them to reconstruct the Lagrangian admitting certain 
solutions. Interestingly, some of the existence conditions do not contradict with each other, 
so the corresponding solutions could coexist in the same Lagrangian. 

As a simplified example, let us construct a model that admits a de Sitter solution, an anti- 
de Sitter solution and a Minkowski solution at the same time. According to the existence 
condition in subsection IVI B[ this can be realized by designing a function f(T) obeying 



2Tf >T - f + T oc T(T + 2Ai)(T + 2A 2 ). 
Assuming the polynomial form of f(T), we build a model 



T 3 + — (Ai + A 2 )T 2 + 20AxA 2 T 



T, 



(83) 



(84) 



where G is a constant with the dimension of mass 2 . When AiA 2 < 0, this model has the 
required multiple solutions indeed. 
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VII. DISCUSSION 



Based on a special frame fl23l) . we found some static solutions with spherical symmetry in 
f(T) gravity theories, where the Maxwell term was taken into consideration. The solutions 
and their existence conditions were summarized in section IVT1 But there are several problems 
unsolved. 

First, the tetrads f lT5|) for FLRW spacetime are widely used in the literature on f(T) 
cosmology. The frame (1231) we employed is the close cousin of (TL31) . However, as we have 
shown in section [Til the reconstruction of Lagrangian is very sensitive to the choice of 
frame, because the torsion scalar T is not locally Lorentz invariant. This poses the pressing 
problem: why are the frames like (Tl3|) so special and what shall we interpret the local Lorentz 



transformations in general f(T) theories? As an up-to-date reference, [40J provides a nice 
answer to this problem for FRLW spacetime, which is further illustrated in [411 ] . 

Second, our ansatz of tetrads simplifies the computation greatly, but it also limits 
the validity of our results. This means we have only found some spherically symmetric static 
solutions for certain special f(T) models. It is still an open problem to get such solutions 
for more general models in other frames. That would involve six more variables related to 
the local Lorentz transformation, leading to equations of motion hard to solve. 

Third, in subsection IV Al we only laid out some particular solutions to T r = 0. Investi- 
gation is still needed to exhaust all of its solutions. 

Fourth, since the existence of black solutions has been confirmed in this paper, the ther- 
modynamics of gravity is to be studied in f(T) theories. 
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